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1.  Introduction 

An  integral  part  of  any  commercial  activity  is  the  periodic  counting 
of  assets.  In  some  larger  establishments  (such  as  a complex  of  warehouses 
and  arsenals)  this  is  a major  undertaking  involving  many  man-hours  of  labor. 
It  is  customary  to  take  such  an  inventory  on  an  annual  basis,  particularly 
for  auditing  purposes.  This  procedure  has  two  major  disadvantages:  (a)  Since 
the  inventory  will  usually  comprise  a very  large  number  of  individual  items 
that  have  to  be  counted,  an  annual  inventory  will  dislocate  normal  procedures 
and  often  require  special  labor;  (b)  if  no  "on-the-spot"  inventories  are 
taken  for  a whole  year,  the  inventory  record  is  liable  to  get  "out-of-date" 
with  a serious  accumulation  of  errors  which  makes,  in  certain  cases,  the 
inventory  records  useless  for  an  up-to-date  picture  on  the  item  counts. 

It  was  with  these  considerations  in  mind  that  M.  R.  Bryson  (1960) 
suggested  a method  of  taking  inventories  piecemeal  on  a sampling  basis.  The 
present  approach  differs  from  Bryson's  method  in  that  it  employs  a sampling 
design  in  which  items  are  inspected  on  a rotating  basis  but  also  uses  the 
time  series  of  inventory  records.  The  present  procedure,  therefore  provides 

(1)  a quality  control  of  the  inventory  records  protecting  the 
establishment  from  having  to  use  inventory  records  that  are 
seriously  in  error,  and 

(2)  a method  of  computing  best  estimates  of  complete  inventories 
for  items  at  any  specified  time  with  computable  statistical 
errors  of  inventory  estimation. 

As  indicated  above,  the  method  is  a combination  of  rotation  sampling 
of  items  and  the  application  of  an  autoregressive  time  series  analysis 
applied  to  the  inventory  records. 


The  method  as  described  in  the  subsequent  sections  applies  to  a situa- 
tion (often  encountered)  where  storage  facilities  are  at  a premium  so  that 
large  incoming  shipments  as  well  as  large  quantities  of  stock  on  hand  are 
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to  be  avoided  if  possible.  Reorder  point  policies  which  allow  for  planned 
large  reorders  and  large  cross  leveling  actions  are  briefly  discussed  in 
Section  8. 

At  the  present  time,  the  process  here  described  could  not  be  applied 
to  actual  inventory  data  since  this  would  have  required  the  adoption  of  our 
technique  as  an  inventory  control  strategy.  However,  it  is  hoped  to  test 
the  procedure  on  simulated  data  in  the  future  before  submitting  the  report 
for  publication. 

2.  Description  of  the  Problem  and  Notation 

Our  notation  will  be  in  accordance  with  Bryson  (1960).  Within  a ware- 
house, or  supply  center,  there  are  different  types  of  items  being  stored. 
Henceforth,  the  term  "item"  will  denote  a type  of  item,  such  as  a bolt  or  a 
more  complex  assembly,  and  the  term  "piece"  will  denote  an  individual  unit 
of  the  type  of  item.  Thus  we  are  interested  in  the  number  of  pieces  of  each 
item  that  are  on  hand  at  a given  time. 

A continuous  record  is  kept  of  all  activities  and  the  resulting  changes 
in  the  amount  of  each  item  on  hand.  Let  denote  this  inventory  record 
for  item  i at  time  t.  The  time  t may  be  measured  in  weeks  or  months  depend- 
ing on  circumstances.  If  a physical  count  is  made  on  item  i at  time  t,  this 
actual  inventory  is  denoted  . 

The  recording  of  activities  for  an  item  is  not  instantaneous.  In  fact, 
due  to  delays  caused  by  paperwork  and  other  administrative  processing,  the 
inventory  record  is  probably  somewhat  out  of  date  compared  with  the  actual 
inventory. 

The  supply  center  is  to  be  "stratified"  so  that  items  of  similar  activ- 
ity and  having  similar  administrative  delays  of  record  keeping  are  placed 
within  the  same  stratum.  We  now  confine  our  discussion  to  activities  within 
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a single  stratum,  noting  that  the  technique  to  be  described  applies  to  all 
strata  on  an  individual  basis. 

3.  Rotational  Design 

Let  us  define  a "short  time  interval"  such  that  the  inventory  records 
are  examined  at  the  beginning  of  each  such  interval.  For  fast  moving  items 
this  interval  may  be  a week.  For  slower  moving  items  it  may  be  a month  but 
should  be  the  same  for  items  in  a stratum. 

The  items  within  a stratum  are  randomly  divided  into  rotational  groups, 
such  that  a group  of  items  is  physically  inventoried  in  a particular  week 
(month)  and  then  is  not  inventoried  again  until  all  such  rotational  groups 
have  been  inventoried.  Therefore,  the  desired  frequency  of  a physical 
inventory  of  each  item  must  be  considered  in  determining  the  number  and  hence 
the  size  of  rotational  groups.  If,  for  example,  it  is  desirable  to  inventory 
each  item  at  least  annually,  there  should  be  no  more  than  52  (12)  rotational 
groups,  perhaps  of  equal  size. 

As  an  illustration  consider  the  rotational  monthly  inventory  plan  for 

twenty-four  items  over  a two  year  period  shown  in  Figure  1. 

We  now  confine  the  scope  of  our  discussion  to  an  individual  item,  noting 

of  course  that  the  discussion  applies  to  all  items  on  an  individual  basis. 

Let  the  term  "period"  denote  the  time  interval  between  successive 

inventories  of  an  item.  Then  period  p is  the  time  interval  which  includes 

the  pth  inventory  of  an  item  and  all  weeks  (months)  preceding,  but  not 

s t 

including,  the  (p  + 1)  inventory  of  the  same  item. 

4.  A Model  for  the  "Normal"  Fluctuation  in  the  Inventory  Record 

The  following  rationale  was  used  to  derive  the  autoregressive  time 
series  model  which  is  described  below.  It  Wl  not  considered  realistic 
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to  assume  that  the  inventory  records  constitute  a stationary  time  series 
for  there  may  well  be  both  upward  and  downward  linear  or  higher  order 
trends  over  the  period  or  (indeed)  seasonal  fluctuations.  However,  follow- 
ing Box  and  Jenkins  (1970)  it  was  assumed  that  a stationary  process  could 
be  reached  if  the  original  time  series  of  inventory  records  would  be  differenced. 
A third  order  difference  was  therefore  assumed  to  constitute  a stationary 
process.  Such  an  assumption  would  allow  for  linear  and  parabolic  like  trends 
over  a period  p. 

To  fix  the  idea,  the  third  difference  in  the  records  for  item  i at  time 
t in  period  p is  given  by 


Zipt  = A3xipt  = (xiPt  ' 3xipt-l  + 3xiPt-2  ' XiPt-3)  • 


(1) 


where  x^  denotes  the  t1-*1  record  of  period  p and  the  third  difference  z^pt 
is  defined  by  (1)  above.  An  assumption  to  be  made  at  this  point  is  that 
under  normal  conditions  the  third  differences  form  a stationary  autoregressive 
process  of  the  first  order.  However,  when  an  inventory  is  taken  on  item  i and 
the  record  is  updated,  there  may  be  a sizable  jump  in  ^ and  hence  the  first 
third  difference  assumed  to  follow  the  stationary  process  is  zip^*  Therefore, 
z^pt  can  be  calculated  only  after  the  4th  weekly  (monthly)  record  is  available 
for  a particular  period. 

Since  our  discussion  is  confined  to  an  individual  item,  for  notational 
simplicity  we  drop  the  item  subscript,  i,  and  define  z to  be  the  third 
difference  at  time  t within  period  p.  And  the  number  of  weeks  (months)  in 
the  p1"*1  period  will  be  denoted  N^. 
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5.  Monitoring 


For  the  purpose  of  monitoring  the  monthly  record,  a time  series  forecast 

of  z ...  is  obtained  from  which  a conditional  confidence  interval  for  x 

pt+1  pt+1 

is  constructed.  When  the  next  record  xpt+^  becomes  available,  it  is  compared 
with  the  confidence  interval  to  determine  whether  or  not  further  action  is 
necessary.  A detailed  discussion  of  this  monitoring  procedure  follows. 

A summary  flowchart  of  the  procedure  to  be  discussed  is  given  in  Figure 
(A.l)  and  Figure  (A. 2)  Appendix  A,  where  the  stages  of  the  procedure  have 
been  numbered  in  a logical  sequence.  In  what  follows  we  shall  frequently 
make  reference  to  this  flowchart. 

To  fit  the  first  order  autoregressive  process 


zpt+l  ^pZpt  + £pt+l  ’ 


an  estimate  of  ^ is  obtained  using  the  method  discussed  by  Box  and  Jenkins 
(1970).  The  computational  formula  for  this  estimate,  denoted  r^(p),  for  a 


completed  period  is 


^(p) 


[tz5(2Pt-rV(zPt“V/(V5)1 

tz  (Zpt-zp)2/(y4)] 


And  for  a period  not  yet  completed,  let  N*  denote  the  number  of  months  for 
which  the  inventory  records  are  available.  Then  substituting  N*  for  N in 

p p 

(3)  gives  the  desired  result.  Note  that  N*  must  be  greater  than  five  in 
order  to  obtain  r^(p). 

A forecast  of  zpt+^  is  then  given  by 


rl(P)  *pt  . t > 6 . 
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The  corresponding  forecast  of  x ...  is  obtained  by  substituting  the  expres- 

pt+i 

sion  for  z from  (1)  into  (4)  above.  This  yields 


x = z _ + 3x  . - 3x  . + x . 
pt+1  pt+1  pt  pt-1  pt-2 

To  obtain  the  conditional  confidence  interval  for  x given  x .,  x , , 

pt+1  pt  pt-i 

x and  x ^ _ we  note  that  the  only  variable  in  (5)  is  the  forecast  of 

pt-/  pt-J 

z that  is  z so  that  the  variance  of  z ...  given  the  preceding  z 

pt+1  pt+1  pt+1  pt 

is  also  the  conditional  variance  of  xpt+^  given  the  preceding  xpt  To 
obtain  the  former  we  use  standard  formulas  for  an  autoregressive  process. 
This  is  given  by 

2 2 '2 

s„  = (1  - rf (P))°y  , 

p i z 

P 

*2  2 

where  a is  the  estimate  of  o , the  variance  of  the  z for  period  p, 
z z pt 

P P 

having  the  computational  form 


(5) 


(6) 


*2 


Nr 


o ~E(z-z)/(N  -4). 

zp  t-4  pt  P P 

Again  we  note  that  N*  > 4 is  substituted  for  when  the  period  is  not 

complete . 

Obtaining  u^  a/2)  ^rom  tlle  Standard  Normal  Tables,  a 100(1  - a)%  con- 
fidence interval  for  Xpt+^  is  given  by 


x ^ ± / s2  u, 


"pt+1  1 ' sp  (l-a/2)  • 


(7) 


We  note  at  this  point  that  during  the  early  weeks  (months)  of  a period, 

the  number  of  third  differences  may  be  insufficient  to  obtain  "reliable" 

2 

estimates  of  <J>  and  o 
P z 


In  fact,  when  months  are  used  as  time  intervals, 


the  number  of  months  in  an  entire  period  may  be  considered  insufficient.  For 
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this  reason  it  is  suggested  that  pooled  estimates  of  the  overall  <j>  and 

be  used  for  the  purpose  of  constructing  confidence  i 'ervals  for  Xpt+^ 

during  the  early  weeks  (months)  of  a period.  If  it  can  be  established 
2 

that  d>  and  a are  constant  for  all  periods,  the  pooled  estimates  might 
P z 

P -2 

be  preferred  over  r^(p)  and  a ^ under  any  circumstances.  However,  a word 

ZP 

of  caution  is  in  order  concerning  the  use  of  the  pooled  estimates  for 

other  than  the  early  weeks  (months)  of  a period,  especially  when  4>  and 
2 

a are  not  constant  over  all  periods.  Bias  may  be  introduced  through  the 
ZP 

pooled  estimates,  which  may  have  such  an  adverse  affect  on  the  forecast 
that  the  whole  purpose  for  pooling  is  lost. 

The  pooled  estimate  of  4>  would  consist  of  a weighted  average  of  the 
r^(p)'s  over  all  periods.  Specifically, 

E(N  -4)r1(p) 

r = 2_ (E 

1 E(N  -4)  ’ U 

P P 

where  r^(p)  for  the  current  period  is  not  available  unless  N*  > 5. 

Thus  using  r^  allows  forecasting  of  Zpt+^  in  the  current  period  as 
soon  as  the  fourth  record  is  available.  Whereas  using  r^(p)  would  require 
that  at  least  the  first  six  records  be  available. 

Letting  r^  replace  r^(p)  in  (4),  the  forecast  of  zpt+^  becomes 


z . ..  = r,z  , t > 4 . 

pt+1  1 pt  — 


The  resulting  forecast  of  xpt+^  is 


Xpt+1  = Zpt+1  + 3xpt  " 3xpt-l  + Xpt-2 
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The  estimate  of  the  conditional  variance  of  x +j  Riven  the  preceding 


x , is  given  by 
pt 


~2  ..  -2'2 
s = (1  - r )o  , 
I z 


where  is  the  pooled  estimate  of  the  variance  of  z over  all  periods  and 
is  given  by 

[E(N  -4)a*  ] 

'2  p P £ 

°z  [E(Np-4)]  ’ 

where  again  N*  replaces  Np  for  the  current  period. 

The  100(1  - a)%  confidence  interval  for  x ..is  now 

pt+i 

Xpt+1  1 /*f  u(l-a/2)  ’ (11) 

where  is  obtained  from  the  Standard  Normal  Tables,  as  before. 

Once  a confidence  interval  for  has  been  constructed,  using  either 

(7)  or  (11)  above,  the  monitoring  test  is  carried  out  by  comparing  the 

observed  x with  the  confidence  interval.  This  can  begin  as  early  as 

the  fifth  month  of  the  current  period.  This  is  step  2 of  Figure  (A.l). 

If  x ,,  falls  outside  the  confidence  interval  it  is  an  indication 
pt+1 

that  either 

(a)  there  is  a sizable  error  in  the  last  record  x 

pt+1 


(b)  an  unexpectedly  large  transaction  has  taken  place  during  the  last 
' interval . 

When  investigating  the  cause  of  the  discrepancy,  (a)  should  be  con- 
sidered first.  If  an  error  is  found  in  the  records,  it  can  be  corrected. 
The  monitoring  test  is  then  repeated  using  the  corrected  records.  If  x . 
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is  still  outside  the  confidence  interval  after  correction,  or  if  no  error 
was  found  in  the  records,  (b)  is  considered.  This  requires  that  an  inventory 
be  taken  to  verify  this  abnormally  large  change  in  the  record.  This  corre- 
sponds to  steps  3 and  4 of  Figure  (A.l). 

When  an  inventory  occurs  due  to  monitoring,  the  item  is  automatically 
assigned  to  the  rotational  group  for  which  the  regular  inventory  is  being 
taken  during  the  present  month.  Thus  an  inventory  for  the  item  is  not  taken 
again  until  its  newly  assigned  rotational  group  is  scheduled  again,  unless 
otherwise  dictated  by  future  monitoring.  This  is  step  5 in  Figure  (A.l). 

6.  The  Estimates  of  Item-Inventories  Under  'Stationary  Conditions' 


It  is  usual  practice  to  use  the  last  inventory-record  of  each 

item  available  at  time  t as  the  best  source  of  information  concerning  the 
inventory  of  item  i in  question.  However,  it  is  reasonable  to  ask  whether 
a better  "update"  of  the  inventory  could  be  estimated  from  the  data.  This 
is  particularly  relevant  if  is  known  to  be  "out  of  date"  owing  to 

unavoidable  delays  in  updating  the  records. 

The  procedure  developed  below  assumes  that  the  third  differences 
of  the  inventory  record  series  x^  follows  the  postulated  stationary 
process.  The  monitoring  of  this  process  as  described  in  Section  3 has  there- 
fore a double  purpose,  namely  (a)  to  invoke  the  taking  of  an  actual  inventory 
y^pt  if  a departure  from  the  stationary  process  is  discovered  and  (b)  the 
prediction  of  y ^ if  the  stationary  process  cannot  be  rejected.  It  is  this 
latter  use  which  we  describe  in  this  section  which  is  step  6 of  Figure  (A.l). 

Figure  2 illustrates  the  method 


’ 
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t+e  t+i 


Figure  2 


In  Figure  2,  the  delay  in  the  inventory  records  is  denoted  9,  where 


0 <_  0 <_  1.  However,  when  weeks  are  used  as  intervals,  it  seems  possible 
that  the  delay  could  exceed  one  week.  In  which  case  the  upper  limit  of 
the  restriction  on  0 may  be  relaxed  with  caution  but  should  not  be  allowed 
to  exceed  two  weeks.  If  in  fact  the  delay  is  more  than  1.5  weeks,  it  may 
be  more  appropriate  to  use  months  as  time  intervals  in  the  future.  Thus 
the  actual  inventory  at  time  t,  y , can  be  written  in  terms  of  0 and  the 
inventory  records  using  either  linear  or  quadratic  interpolation.  We  now 


consider  the  estimation  of  0 and  the  resulting  forecast  of  ypt  for  both 


cases . 
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The  linear  interpolate  is  of  course  the  simpler  of  the  two.  It  is 


given  by 


ypt  = (1  - 0)xpt  + 0xpt+l  ’ 0 - 6 - 1 • 


And  for  the  purpose  of  obtaining  an  estimate  of  0,  this  expression  can  be 
rewritten  as 

(ypt  - xpt)=  e(xpt+i  - V * 0<e<i.  (13) 

It  should  be  recalled  at  this  point  that  the  items  have  been  stratified 
such  that  all  items  within  a stratum  have  values  of  0 which  can  be  assumed 
essentially  equal.  Therefore,  it  is  possible  to  obtain  an  estimate  of  0, 
say  0L>  using  least  squares,  by  fitting  (13)  to  the  available  data  Within 
a stratum. 

Returning  now  to  the  triple  subscripting,  the  data  will  include  y ^ 

and  the  corresponding  records,  x^^  and  xip2>  for  all  i and  p,  except  where 

the  inventory  was  taken  due  to  monitoring.  Note  that  since  a forecast  of 

x.  cannot  be  obtained  until  t = 4,  the  record  x.  „ must  be  available 
ipt+1  ip2 

before  y.  , can  be  included  in  the  data, 
ipl 

After  computing  0 , the  restriction  0 < 0 < 1 is  imposed.  When  months 

L L 

are  being  used  as  intervals,  if  0 < 0,  use  0 = 0;  or  if  0T  > 1,  use  0 = 1. 

As  mentioned  previously,  the  upper  limit  may  be  somewhat  relaxed  when  weeks 
are  used  as  intervals.  Thus,  in  the  case  where  weeks  are  used,  if  0,  < 0 use 

li 

0^  ” 0,  but  if  0j  > 1,  0L  may  either  be  used  with  caution  or  set  equal  to 
unity.  This  will  obviously  give  the  restricted  least  squares  estimate  of  0. 

The  forecast  of  the  current  inventory  is  then  obtained  by  substituting 
0^  for  0 and  x^pt+j  f°r  xipt+i  *-n  (12).  This  gives 
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yipt  = (1  - eL)xipt  + Vipt+l 


where  x.  , . Is  the  forecast  of  x.  ,,  from  either  (5)  or  (10). 
ipt+1  ipt+1 


When  the  time  interval  is  a month,  it  may  be  more  appropriate  to  use 


a quadratic  interpolation  formula.  The  quadratic  interpolate  is  given  by 


ylpt  = %0(1  - + (1  - 02)xipt  + %6(1  + e)*ipt+1 


0 < 0 < 1 


The  estimate  of  0,  say  0^,  is  obtained  from  the  available  data  in  a stratum 


as  follows. 


Define 


Q(9)  * l ' ty(i,P,l)  - %9(1  - 6)x(i,p-l,N  ) 


(1  - 02)x. v - %0(1  + 0)x,. 


(ipl) 


(ip2) 


Q (0 ) is  tabulated  as  a function  of  0 (say  at  interval  A0  = .05),  for 


0 < 0 < 1.2.  The  absolute  minimum  is  then  determined  numerically,  using 


Mathematical  Programming  by  Scanning.  The  value  of  0 corresponding  to  this 


absolute  minimum  is  then  the  desired  estimate,  0^. 


The  forecast  of  the  current  inventory  is  then  obtained  by  substituting 


0p  for  0 and  xipt+^  f°r  xipt+i  in  0-5).  This  yields 


yipt  * %0Q(1  ‘ 0Q)xipt-l  + (1  “ 9Q)xipt  + %9Q(1  + 0Q)xiPt+l 


where  x.  is  the  forecast  of  x.  from  either  (5)  or  (10). 
ipt+1  ipt+1 


Since  new  data  will  become  available  each  week  (month),  the  estimate 


of  0,  either  0T  or  should  be  updated  frequently. 
L W 
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7.  A Quality  Control  Chart  Monitoring  Regular  Inventories  Against  Inventory 


Records 


When  the  p inventory  is  taken  on  item  i,  a new  period  p begins  with 
as  the  updated  value  of  the  inventory  record.  (Step  1,  Figure  (A. 2).) 
A discrepancy  between  y and  xipi  *s  expected  due  to  the  delay  in  the 

records.  However,  if  this  discrepancy  is  unusually  large,  it  may  be  an 
indication  that  either 

(a)  errors  exist  in  the  records  which  may  have  accumulated  during  the 
previous  period 

or  in  the  case  where  y is  extremely  smaller  than 

(b)  pilfering  or  other  unexplained  losses  have  occurred. 

In  order  to  monitor  such  discrepancies,  a control  chart  technique  can 


be  implemented.  Define 


d.  = y.  , - x , 
ip  ipl  ipl 


where  only  the  regular  inventories  (not  resulting  from  the  monitoring 
discussed  in  Section  5)  are  included.  Recalling  again  that  the  items  have 
been  stratified  according  to  activity  and  delay,  an  estimate  of  the  variance 
of  the  discrepancies  is  given  by 

s.  = Eld  /n  , 
d iP  fI> 

where  n is  the  total  number  of  d.  's  available  in  the  stratum.  100(1  - a)% 

iP 

control  limits  are  then  given  by 


1 “(l-a/2)  ' 


where  u 


(l-c/2) 


(Step  2,  Figure  (A. 2).) 


is  from  the  Standard  Normal  Tables. 
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If  d . falls  outside  the  control  limits,  (a)  is  considered  first, 

ip 

If  an  error  is  found  in  the  record,  it  can  be  corrected.  If  after 

correcting  the  records  d^  f^lls  below  the  lower  control  limit,  or  if  no 

error  is  found  and  d.  is  below  the  lower  limit  (b)  is  considered.  In 

ip 

this  case  appropriate  action  should  be  taken  to  prevent  further  unex- 
plained losses.  (Steps  3,  4,  and  5 of  Figure  (A. 2).) 
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APPENDIX  A 

Figure  (A. 1) 

Summary  Flowchart 


cont . 
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Figure  (A. 2) 


I 


*‘v  ^ * 


-17- 


8.  Pre-Planned  Re-Orders  and  Cross-Leveling  Actions 

As  mentioned  in  Section  1,  the  situation  covered  by  the  above  proce- 
dure is  one  (often  encountered)  where  storage  facilities  are  at  a premium 
so  that  large  incoming  shipments  as  well  as  large  quantities  of  stock-on- 
hand  are  to  be  avoided  if  possible.  It  is  realized  that  this  is  not  always 
feasible  or  practical  and  the  following  modification  of  the  procedure  may 
then  be  implemented: 

In  situations  where  large  changes  in  the  inventory  are  planned  as  part 
of  the  normal  reorder-point  policy,  management  may  distinguish  between 
(a)  "expected  or  planned  sizeable  changes  in  the  inventory"  and  (b)  "unex- 
pected sizeable  changes  in  the  inventory." 

In  the  case  of  (a),  after  checking  the  correctness  of  the  inventory 
records,  this  may  be  accepted  as  correct  without  insisting  on  an  on-the- 
spot  inventory.  However,  it  would  be  necessary  to  restart  the  time  series 
analysis  accepting  the  new  record  as  a starting  point. 

In  the  case  of  (b) , it  is  suggested  that  the  present  procedure  be 
adopted  since  unexpected  changes  in  the  inventory  (particularly  inventories 
below  a safety  level)  should  be  subjected  to  an  on-the-spot  check. 
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